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The superconductor-insulator transition of ultrathin films 
of bismuth, grown on liquid helium cooled substrates, has 
been studied. The transition was tuned by changing both 
film thickness and perpendicular magnetic field. Assuming 
that the transition is controlled by a T — critical point, a 
finite size scaling analysis was carried out to determine the 
correlation length exponent v and the dynamical critical ex- 
ponent z. The phase diagram and the critical resistance have 
been studied as a function of film thickness and magnetic field. 
The results are discussed in terms of bosonic models of the 
superconductor-insulator transition, as well as the percolation 
models which predict finite dissipation at T — 0. 



I. INTRODUCTION 

After about two decades of 

research, the superconductor-insulator (SI) transition in 
disordered films of metals remains a controversial sub- 
ject, mainly due to contradictory results in both the- 
oretical and experimental studies. This work aims to 
improve the understanding of this phenomenon, which 
might also be relevant for high-T c superconductors and 
possibly connected to novel metal-insulator transitions in 
2D electron systems. 

The superconductor-insulator transition in ultrathin 
films of metals is believed to occur at the absolute zero of 
temperature when the quantum ground state of the sys- 
tem is changed by tuning disorder, film thickness, carrier 
concentration or magnetic field. Unlike finite tempera- 
ture phase transitions in which thermal fluctuations are 
crucial, T = phase transitions are driven purely by 
quantum fluctuations. At finite temperatures, an un- 
derlying quantum phase transition manifests itself in the 
scaling behavior of the resistance with the appropriate 
tuning parameter and the temperature, along with the 
coherence length and dynamical critical exponents, v and 
z respectively Q . Various models of the superconductor- 
insulator transition in disordered films can be roughly 
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divided in two groups: those in which the superconduc- 
tivity is destroyed by fluctuations of the amplitude of the 
order parameter, and those which focus only on the phase 
fluctuations. 

If superconductivity is destroyed only by phase fluctu- 
ations, then Cooper pairs persist on the insulating side 
of the transition and the transition may be described by 
a model of interacting bosons in the presence of disor- 
der. Based on this assumption, Fisher and co-workers 
p| suggested a scaling theory and a phase diagram for 
a two-dimensional system as a function of temperature, 
disorder, and magnetic field The superconducting 

phase is considered to be a condensate of Cooper pairs 
with localized vortices, and the insulating phase is a con- 
densate of vortices with localized Cooper pairs. At the 
transition, both vortices and Cooper pairs are mobile as 
they exchange their roles, which leads to a finite resis- 
tance. Some important predictions of the model are the 
universal value of this critical resistance and specific val- 
ues of the critical exponents v and z. 

This so-called "dirty boson" problem has been ex- 
tensively studied using quantum Monte Carlo simula- 
tions |p|-|l2| , real-space renormalization group techniques 
strong-coupling expansions pq| and in other ways 
Finite temperature behavior in the vicinity of 
a quantum critical point was also studied analytically 
pl"| , p2[ . A transition from a superfluid to a Mott insu- 
lator was found in the pure case, and to a Bose glass 
insulator in the presence of disorder, but there is still 
considerable disagreement as to the universality class of 
the transition, as well as the value of the critical resis- 
tance. 

An alternative picture of interacting electrons (2^-|||] 
proposes a different mechanism: the density of states and 
the Cooper pairing are suppressed on the insulating side 
of the superconductor-insulator transition due to an en- 
hanced Coulomb interaction. The SI transition occurs 
as a consequence of fluctuations in the amplitude, rather 
than the phase of the order parameter. In other words, 
Cooper pairs break up into single electrons at the tran- 
sition. Therefore the superconducting gap would also 
vanish at the transition. 

The model of interacting electrons has also been stud- 
ied numerically. Quantum Monte Carlo simulations of 
an attractive fermion Hubbard model with on-site inter- 
actions p6f yielded a direct superconductor-to-insulator 
transition in two dimensions without an intervening 
metallic phase. The critical resistance was found to de- 
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pend on the strength of the attractive interaction, as 
a function of which, a crossover from a fermionic to a 
bosonic regime occurs. The results of this theory quali- 
tatively resemble the experimental data. A recent calcu- 
lation of the effect of disorder on the gap in the density 
of states, using a similar model P7I , showed that the ex- 
istence of a gap on the insulating side of the transition 
depends on the coupling strength, allowing for a Fermi 
insulator at weak and a Bose insulator at strong coupling. 

Experimentally, the destruction of superconductivity 
by disorder has been studied in films of MoGe p8| , p9[ , 
InO x |0|-|3 and Bi, Pb, Ga, Al ff-|7| among oth- 
ers. Evidence was found of T c going to zero with in- 
creasing disorder |2£| implying the destruction of Cooper 
pairs at the transition. Tunneling experiments also seem 
to support the fermionic picture. Valles et al. found 
that the superconducting gap and the mean field transi- 
tion temperature are both suppressed as disorder is in- 
creased, and that the gap vanishes on the insulating side 
of the superconductor- insulator transition |58|. Hsu et 
al. carried out tunneling studies of the superconductor- 
insulator transition in PbBi/Ge films, and found a large 
number of quasiparticle states near the Fermi energy [ p9[ . 
They estimated the average number of Cooper pairs in 
a coherence volume to be on the order of one at the 
superconductor-insulator transition. This result, in com- 
bination with the disappearance of the energy gap, was 
interpreted as evidence of the superconductor-insulator 
transition being driven by fluctuations in the amplitude 
of the order parameter. Alternatively, it is possible for 
the superconducting energy gap to be reduced or the 
tunneling density of states to be broadened as a con- 
sequence of phase fluctuations |4(j]. Thus, the absence of 
the gap in these tunneling studies does not necessarily 
mean that Cooper pairs are absent on the insulating side 
of the superconductor-insulator transition, but it may 
imply that a full picture might have to include fermionic 
degrees of freedom. 

Evidence of the importance of the bosonic picture can 
be found in the work of Paalanen et al. |u| . These work- 
ers studied the magnetoresistance and the Hall effect in 
amorphous InOa; films and observed two distinct transi- 
tions: one at a critical field B xx where the longitudinal 
resistance diverges and the system presumably undergoes 
a transition from a superconducting phase to a Bose glass 
insulator with localized Cooper pairs, and the other at a 
higher field B° y , where the transverse resistance diverges 
and the Cooper pairs of the Bose glass insulator presum- 
ably unbind. The transition in the transverse resistance 
occurred at the same magnetic field where the longitudi- 
nal resistance showed a maximum. Since a Bose insulator 
might be expected to have a higher resistance than an in- 
sulator with localized single electrons, and from the dis- 
order dependence of B xx /B xy , this was interpreted as ev- 
idence of the bosonic nature of the insulating state close 
to the superconductor-insulator transition. Similar be- 
havior was observed by other groups Q . Magnetoresis- 
tance studies of amorphous InO x films by Gantmakhcr 



et al. |32j also seem to support the bosonic picture. Fur- 
thermore, a linear component of the magnetoresistance 
observed in the insulating regime in amorphous Bi films 
can be interpreted as a signature of vortex motion [fi~l| . 

In the context of the scaling behavior, the thickness 
tuned transition of ultrathin films of amorphous Bi has 
been studied in zero magnetic field |3(| . A scaling anal- 
ysis of the magnetic field tuned SI transition has been 
carried out for thin films of InO x @ and MoGe. |j| 
All of these investigations found v ~ 1.3 and z « 1, 
consistent with the theoretical predictions of the boson 
Hubbard model. 

Yet another interpretation of the experimental data 
has recently been proposed by Shimshoni et al. and 
expanded upon by Mason and Kapitulnik ^3|. In this 
picture, a film contains both insulating and supercon- 
ducting puddles, and transport is dominated by tunnel- 
ing or activated hopping between them. The SI transition 
then occurs as a consequence of the percolation of one 
phase or the other. Since the correlation length exponent 
in 2D classical percolation is 4/3, this is consistent with 
f«1.3 observed in most experiments. This model also 
predicts a saturation of the resistance at very low temper- 
atures, which seems to be supported by the experimental 
data of Ephron et al. [Q, and Yazdani and Kapitul- 
nik p9| . Similar effects have been observed in the much 
earlier work of Wang et al. [ |i"5| on underdoped high-T c 
(cuprate) films. These ideas may be relevant to similar 
features of the results of Kravchenko et al. on two dimen- 
sional electron gas systems [Q. In all studies in which 
there is flattening in R(T) at low temperatures, one must 
be concerned with the possibility of electrical noise being 
the source of the effect. Also in multi-component ma- 
terials such as MoGe and underdoped cuprates there is 
always a possibility of second phases affecting the out- 
come. Furthermore, it has recently been proposed that 
the flattening in R(T) at low temperatures may be a sig- 
nature of Bose metal, a phase in which the Cooper pairs 
are mobile but do not condense (47| . 

The quantitative results of the study of the magnetic 
field tuned superconductor-insulator transition presented 
here for disordered metal systems are in serious dis- 
agreement with previous measurements of this transition, 
adding yet another puzzle to this problem, and calling for 
a re-examination of existing models. The thickness-tuned 
transition has also been studied in a nonzero magnetic 
field. This allows for the construction of a phase dia- 
gram and a direct comparison of the two different ways 
of tuning the SI transition, by varying thickness or mag- 
netic field. 

This paper is organized as follows: the finite-size scal- 
ing procedures used to determine the critical exponents 
are described in Section II. Experimental details are given 
in Section III. Section IV focuses on the magnetic field- 
tuned transition, while the analysis of the thickness- 
tuned transition in finite magnetic field, which has not 
been studied before, is presented in Section V. In Sec- 
tion VI, the phase diagram as a function of thickness 
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and magnetic field is presented. The critical resistance 
and its apparent non-universality are discussed in Section 
VII. The results and their implications are summarized 
and further discussed in Section VIII. A brief account of 
a portion of this work has been previously reported |M8| . 



II. SCALING PROCEDURES 

The scale of fluctuations on either side of a quan- 
tum phase transition is set by a diverging correlation 
length £ oc 5~ v and a vanishing characteristic frequency 
51 oc £ -z . Here S is the deviation from the critical point 
S = \K—K C \, where K is the control or tuning parameter, 
which drives the system through the transition (i.e. dis- 
order, thickness, magnetic field, etc.), K c is the critical 
value of K at the transition, v is the correlation length 
exponent and z is the dynamical critical exponent. The 
exponents v and z determine the universality class of the 
transition. They may not depend on the microscopic de- 
tails of the physics of the system under study, but on its 
dimensionality, the symmetry group of its Hamiltonian 
and the range of interactions. 

The resistance of a two-dimensional system in the 
quantum critical regime follows the scaling relation jl],[3| : 



R(S,T) = R C fiST- 1 ^ 2 ) 



(1) 



Here S = \d — d c \ in the case of the thickness-tuned tran- 
sition and S = \B — B c \ in the case of the magnetic 
field-tuned transition. R c is the critical resistance at 
6 = 0, and f(x) is a universal scaling function such that 
/(0) 1- 

The first step in the analysis of the experimental data 
is to determine the critical value of the tuning parameter 
and plot the resistance as a function of 5. The <5-axis is 
then re-scaled by a factor t: 



R(S,t)=R c f(6t) 



(2) 



where the parameter t{T) is determined at each tem- 
perature by performing a numerical minimization which 
yields the best collapse of the data. If the resistance re- 
ally follows the scaling law (Eq. |l|), it is obvious that t(T) 
has to be a power law in temperature, i(T) = T~ l ' uz . 
The exponent product vz is then found by plotting t(T) 
as a function of T on a log-log scale, and determining the 
slope which is then equal to — 1/vz. 



Similarly, at a constant temperature 29 



R(S, E) = R c f(6E 



-l/v(z+l) 



(3) 



where E is the electric field across the sample. This 
time, the i5-axis is re-scaled by a field-dependent fac- 
tor, t(E), which should be a power law in electric field, 
t(E) = £ , ~ 1 / l/ ( z + 1 ) , and the exponent v(z+l) can then be 
determined from the field dependence of the parameter 
t(E). 



The main advantage of this scaling procedure is that it 
requires neither prior knowledge of the critical exponents, 
nor the temperature and thickness dependence of the re- 
sistance. The critical exponents are determined empiri- 
cally from the data, with the critical exponent product as 
the only adjustable parameter, while the critical value of 
the tuning parameter is determined independently. The 
temperature scaling determines the product vz, while the 
electric field scaling determines v(z + 1). Combining the 
two results, the correlation length exponent v and the 
dynamical exponent z can be determined separately. 

An alternative way to determine these critical expo- 
nent products is to evaluate a derivative of the resistance 
with respect to K at its critical value K c [M : 



{dR/dK) Kc oc R C T 



-1/vz 



/(0) 



(4) 



where K = d at the thickness-tuned transition and K = 
B at the magnetic field-tuned transition, and /'(0) is a 
constant. Plotting (dR/ 8K)k c as a function of T _1 on 
a log-log scale should yield a straight line, with a slope 
equal to 1/vz. The same method can be applied to the 
electric field scaling to determine l/v(z + 1), and then v 
and z can be calculated from the results. 

In the work described below, both scaling procedures 
were used to obtain the critical exponents, in order to 
check their consistency. The exponents obtained using 
two different methods were found to be the same, within 
the experimental uncertainty. 



III. EXPERIMENTAL METHODS 

Ultrathin Bi films were evaporated on top of a 10A 
thick layer of amorphous Ge, which was pre-deposited 
onto either SrTiO^ or glazed alumina substrates. The 
substrate temperature was kept well below 20K during 
all depositions and all the films were grown in situ un- 
der UHV conditions (~ 1CP 10 Torr). The film thickness 
was gradually increased through successive depositions 
in increments of 0.1 — 0.2 A Resistance measurements 
were carried out between the depositions using a stan- 
dard DC four-probe technique, with currents up to 50 
nA. A detailed temperature dependence of the resistance 
in zero field and in magnetic field was recorded at each 
film thickness in the temperature range between 0.14K 
and 15K, where the lowest temperatures were achieved 
using a dilution refrigerator. As the film thickness in- 
creased from 7A to 15A, the temperature dependence of 
the resistance of the system changed from insulator-like 
to superconductor-like at low temperatures, with no sign 
of reentrant behavior typically observed in granular films 
|p5]| . The films that were superconducting in zero field 
were driven insulating by applying a magnetic field of 
up to 12 kG perpendicular to the plane of the sample 
using a superconducting split-coil magnet. The scaling 
procedures described above were applied to the magnetic 
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field-tuned transition, as well as to the thickness-tuned 
transition in both zero field and in a fixed magnetic field. 



IV. MAGNETIC FIELD-TUNED SI TRANSITION 

The resistance as a function of temperature for seven 
films with varying degrees of disorder was studied in mag- 
netic fields up to 12 kG applied perpendicular to the 
plane of the sample. A typical temperature dependence 
of the resistance as the magnetic field changes is shown 
in Fig. |l|. In zero field, the resistance decreases with 
decreasing temperature suggesting the existence of su- 
perconducting fluctuations. A magnetic field destroys 
this downward curvature, and at some critical magnetic 
field, B c , the resistance is independent of temperature. 
In magnetic fields higher than B c the film is insulating, 
with OR/dT < 0. Figure || shows the resistance as a 
function of magnetic field for different temperatures. 

If the sheet resistance is normalized by the value of 
the critical resistance at each thickness, R/R c (d), then 
all the data can be collapsed onto a single curve. The 
collapse of the normalized resistance data as a function 
of 5t for five samples is shown in Fig. |^. The critical 
exponent product vz, determined from the temperature 
dependence of the parameter t (inset of FigJ|. ) , is found 
to be vz = 0.7 ± 0.2, apparently independent of the film 
thickness. The same exponent products were obtained 
using the alternative method of plotting (dR/ dB)s a vs. 
T~ l on a log-log plot and determining the slope which is 
equal to 1/vz, as shown on Fig. ||. 

Electric field scaling was also carried out for one of the 
samples. Unfortunately, there was not enough data avail- 
able for the insulating side of the transition to carry out 
a complete analysis, but the data on the superconducting 
side was sufficient to obtain the value of the critical ex- 
ponent product v(z + 1). The magnetic field dependence 
of the sheet resistance for different values of electric field 
applied across the sample is shown on Fig. §. The resis- 
tance data were then plotted as a function of [B — B c ), 
and re-scaled by a parameter t{E) to obtain the best 
collapse of the data, shown in Fig. [6l For the electric 
field dependence of the parameter t(E), shown in the in- 
set of Fig. the best power law fit was obtained for 
v(z + 4. Combining this result with the result of the 
temperature scaling, it follows that z « 1 and v fss 0.7 for 
the magnetic field tuned superconductor-insulator tran- 
sition. 

In contrast with our findings, previous studies of thin 
films of amorphous InO x |KJ and MoGe |^9| both 
showed v m 1.3 and z = 1 for the magnetic field tuned 
superconductor-insulator transition. Our surprising re- 
sult is also in obvious disagreement with the prediction 
of the scaling theory (from which v > 1 for a dis- 

ordered system), as well as with the percolation-based 
models |4^] (from which v w 1.3 would be expected). 



V. THICKNESS-TUNED SI TRANSITION 

For very thin films, the resistance increases exponen- 
tially with decreasing temperature, while for the thicker 
films the resistance goes to zero as the films become su- 
perconducting. At the critical thickness d c , the resistance 
is temperature independent, and the system is expected 
to stay metallic down to T = 0. 

Using the same methods described above, the criti- 
cal exponent product vz was determined to be 1.2 ± 0.2 
when the superconductor-insulator transition was tuned 
by changing the film thickness in zero magnetic field [Q . 
A similar scaling behavior has been found in ultrathin 
films of Bi by Liu et al. [Q, with the critical exponent 
product vz = 2.8 on the insulating side and vz = 1.4 
on the superconducting side of the transition. The fact 
that vz was found to be different on the two sides of 
the transition raises the question of whether the mea- 
surements really probed the quantum critical regime. It 
is likely that the scaling was carried out too deep into 
the insulating phase, forcing the scaling form (Eq. [I].) 
on films which were in a fundamentally different insulat- 
ing regime. Such films should not be expected to scale 
together with the superconducting films, hence the dis- 
crepancy on the insulating side of the transition. In the 
present work, the measurements were carried out at lower 
temperatures than previously studied and with more de- 
tail in the range of thicknesses close to the transition. 
Both sides of the transition scaled with vz rj1.2, which 
is close to the value obtained by Liu et al. on the super- 
conducting side of the transition. This result is also in 
very good agreement with the predictions of the scaling 
theory || , renormalization group calculations 0{b||l9| , 
and Monte Carlo simulations 0,^|J|]. 

All previous experiments which studied the thick- 
ness or disorder tuned superconductor-insulator transi- 
tion were carried out in zero magnetic field. An applied 
magnetic field is generally expected to change the univer- 
sality class of the transition since it breaks time reversal 
symmetry. One would therefore expect the critical expo- 
nent product vz to be different in the presence of a finite 
magnetic field. Furthermore, the thickness-tuned transi- 
tion in a finite magnetic field might be expected to be 
in the same universality class as the magnetic field-tuned 
transition at fixed thickness. 

The thickness-tuned superconductor-insulator transi- 
tion in a finite magnetic field was probed by sorting the 
magnetoresistance data which were carefully taken as a 
function of temperature and magnetic field for each film. 
A detailed scaling analysis was carried out at fixed mag- 
netic fields of: 0.5 kG, 1 kG, 2 kG, 3 kG, 4.5 kG and 7 
kG for one set of films, and 12 kG for a different set of 
films. For each value of the magnetic field, the resistance 
was plotted as a function of the film thickness at different 
temperatures, ranging from 0.14 K to 0.5 K, in order to 
determine the critical thickness at that field. If the sheet 
resistance is normalized by the critical value at each field 
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R/R C (B), then the normalized resistance data as a func- 
tion of the scaling variable for all temperatures and all 
values of the magnetic field collapsed onto a single curve, 
as shown on Fig. |t]. The critical exponent product de- 
termined from the parameter t(T) (Inset of Fig. ^) was 
found to be vz = 1.4 ± 0.2, apparently independent of 
the magnetic field. Once again, the alternative scaling 
procedure yielded very similar results, as shown in Fig. 

I 

This value of the product vz is a factor of two larger 
than that obtained for the magnetic field-tuned transi- 
tion. It is, however, very close to that obtained from the 
analysis of the zero-field transition carried out using data 
from the same set of films, which was vz — 1.2 ± 0.2. 
Given the experimental uncertainties, it is hard to say 
whether this difference in value of the exponent products 
reflects a difference between the universality classes of the 
thickness driven transitions in zero and finite magnetic 
field. These exponent products are close to those found in 
Monte Carlo simulations of the (2+l)-dimensional clas- 
sical XY model with disorder by Cha and Girvin |?| . 



VI. THE PHASE DIAGRAM 

Combining the data obtained from the thickness-tuned 
transitions in a fixed magnetic field and the field-tuned 
transitions at the fixed thickness, one can construct a 
phase diagram with thickness and magnetic field as in- 
dependent variables. This is shown in Fig.|| The films 
characterized by parameters which lie above the phase 
boundary are "insulating" (dR/dT < at finite tem- 
peratures), and the ones below it are "superconducting" 
(dR/dT > at finite temperatures). The phase bound- 
ary is a power law: 



B c oe \d — d c 



(5) 



The best fit to the data yields x = 0.7. Near the critical 
thickness for the zero field transition, a simple dimen- 
sionality argument 0j suggests that the critical magnetic 
field should scale as: 



B r OC 



$0 



(6) 



where $o is the flux quantum. Since the correlation 
length is £ cx \d — d c \ , one might expect the critical 
field to be: 



B c oc \d — d c 



(7) 



According to the phase boundary obtained in this exper- 
iment (see Eq. g[), this would mean that v = 0.35, a 
value not consistent with the results of the scaling anal- 
ysis carried out on the same films. It also does not agree 
with v = 1.3 obtained by Refs. @ and @. There is no 
obvious physical reason for such a small value of v and 



implied large values of z, so this discrepancy is a mys- 
tery at this time. It is possible that the simple argument 
expressed in Eqs. [6] and |?] is too naive. 

Another surprising feature of the experimental results 
is that the critical exponent product vz evidently de- 
pends on whether the phase boundary is crossed ver- 
tically (changing the thickness at a constant magnetic 
field), in which case vz ss 1.4, or horizontally (changing 
the magnetic field at a fixed thickness), in which case 
vbzb ~ 0.7. One might expect the critical exponents to 
not depend on the direction in which the boundary is 
crossed. If, however, the actual tuning parameter were 
not film thickness, but some other physical parameter 
which was a function of thickness, a factor of two in 
the critical exponent product determined from an analy- 
sis using thickness rather than the "correct" control pa- 
rameter might result. The "correct" control parameter 
might be some measure of disorder, electron screening, 
damping, or Cooper pair density. The detailed functional 
form of the thickness dependence of these parameters for 
qucnch-condensed films is not known. 

Another possibility is that there are actually two phase 
boundaries, separating three different regimes, so that 
each exponent belongs to a different phase boundary. 
There has been some indication of a vortex liquid phase 
in between the superconducting (vortex glass) phase and 
the insulating (Bose glass) phase [ j50|]5~i| ] . Since there only 
appears to be one phase boundary, that is probably not 
the case. It is possible, however, that the two boundaries 
could be indistinguishable over the range of parameters 
explored in these studies, but would become apparent at 
higher fields, greater film thicknesses, or lower tempera- 
tures. These matters need to be investigated further. 



VII. THE CRITICAL RESISTANCE 

The critical resistance for the field-driven transition, 
contrary to the predictions of the dirty boson models, 
does not seem to be universal. Figure |l0| shows that 
R c decreases as the critical field increases, roughly in 
a linear fashion. Since thicker films have lower normal 
state resistances and higher critical fields, this also means 
that R c decreases with increasing thickness and decreas- 
ing normal state resistance. Very similar behavior was 
observed by Yazdani and Kapitulnik J29| . In order to ex- 
plain the non-universal behavior of the critical resistance, 
these authors proposed a two-channel conduction model, 
in which the conductance due to the electron (fermion) 
channel adds to the conductance due to the boson chan- 
nel. When the unpaired electrons are strongly localized, 
the conduction is mostly due to bosons, and the resis- 
tance is close to Rq = h/Ae 2 , as predicted by the boson 
Hubbard model. In the opposite limit, unpaired electrons 
contribute significantly to the conduction at the transi- 
tion. Films with lower normal state resistances would 
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then have lower critical resistances due to the larger frac- 
tion of normal electrons. The critical resistances in our 
experiment, however, are all greater than Rq and their 
values could only be explained this way if the quantum 
resistance due to pairs was itself greater than Rq . 

The conductance due to the electronic channel in a 
magnetic field might also depend on the strength of the 
spin-orbit interactions, which is another difference be- 
tween our samples and those of Refs. Q and p9[ . The 
strength of the spin-orbit interactions is typically pro- 
portional to Z 4 , where Z is the atomic number. Since 
Bi is a heavy metal, spin-orbit interactions are stronger 
than in the lighter InO x and MoGe. ft is known that 
in the weakly localized systems with strong spin-orbit 
interactions the magnetoresistance is positive, while it 
is negative otherwise 52 5^]. If weakly localized un- 
paired electrons really contributed significantly to the 
conduction at the magnetic field-tuned superconductor- 
insulator transition at the experimentally accessible finite 
temperatures, the contribution to the magnetoresistance 
due to localization effects could have a positive or a neg- 
ative sign, depending on the strength of the spin-orbit 
interactions. This would make R c bigger in the case of Bi 
films, and smaller in the case of InO x and MoGe, consis- 
tent with experimental observations. There is, however, 
a striking similarity in the magnetic field and normal 
state resistance dependence of the critical resistance of 
the Bi films and MoGe films of Ref. J2^]: even though 
their critical resistances fall on the opposite sides of Rq, 
they both decrease with magnetic field roughly linearly, 
with almost the same slope. 

Strictly speaking, the critical resistance is predicted to 
be universal only at T — 0, while the finite temperature 
corrections are expected to be scaled by the with the 
Kosterlitz-Thouless transition temperature, T c ||: 

R C (B C ,T) = R* C +0(J) 2 (8) 

where R* is the universal resistance at T = 0, and R c is 
the critical resistance at some finite temperature as mea- 
sured in the experiments. A closer look at the crossing 
plots such as that of Fig. ||, reveals that the critical resis- 
tance is indeed slightly temperature dependent. A con- 
siderable amount of noise over the accessible temperature 
range made it hard to compare this temperature depen- 
dence with Eq. but qualitative behavior is shown on 
Fig. [O. Normal state resistances of the MoGe films |2^] 
are a factor of 3-10 lower than the Bi films considered 
here, which means that our samples are probing a differ- 
ent part of the phase diagram (normal state resistances 
are inversely proportional to the film thickness in our ex- 
periment), and the finite temperature corrections might 
be more important in one case then the other. Indeed, 
somewhat higher critical resistances were found in InO x 
films if the temperature dependence of R c is taken into 
account [p3[ . 

A recent analytical calculation of the critical resistance 
of a two dimensional system at finite temperatures in the 



dirty boson model including Coulomb interactions |2^] 
yielded a critical resistance of kIARq. The author sug- 
gested that the next order correction would bring the 
result closer to Rq. This result is in excellent agreement 
with the critical resistance found in the present measure- 
ments, which was 1.1 — 1.2Rq. Monte Carlo simulations 
of the (2 + l)-dimensional XY model without disorder Q 
also find the critical resistance to be R c — 7.7fc£l, again 
very close to the value found in this work. 



VIII. DISCUSSION 

A lot of attention has been focused recently on the 

|43Jil|,&|6 



effects of dissipation on SI transitions 42 



Within the picture proposed by Shimshoni et al. [42|, 
the transition between the superconducting and the in- 
sulating state is of a percolative nature. On the insu- 
lating side of the transition, electrical transport occurs 
through activation or tunneling of Cooper pairs between 
superconducting domains. Likewise, on the supercon- 
ducting side, vortices tunnel from one insulating domain 
to another. Using incoherent Boltzmann transport the- 
ory, Shimshoni et al. derive resistivity laws in differ- 
ent temperature regimes and predict finite dissipation at 
T = for all values of the magnetic field. Their re- 
sults seem to be supported by measurements on several 
different systems: thin films E3J29], 2D Josephson junc- 
tion arrays ||, Si MOSFETsHgfand QH systems |§, 
where a saturation of the resistance at low temperatures 
is observed and attributed to dissipation effects. The 
percolative nature of the transition can explain the value 
of v « 1.3 found in most of the field-tuned experiments 
on thin films ]30| , p9[ |, as well as the apparent symmetry 
between insulating and the conducting phase observed in 
other experiments [ ^7||4^j5^ , |59|l . 

In contrast with the above mentioned results, we do 
not observe any saturation in the temperature depen- 
dence of the resistance as the temperature decreases, or 
in other words, SR/ST is non-zero down to the lowest 
temperatures, which were above 0.1K. Of course inves- 
tigation down to even lower temperatures might lead to 
a different conclusion. However a satisfactory fit to the 
resistivity laws predicted by Shimshoni et al. |4^] could 
not be obtained. 

Mason and Kapitulnik recently proposed a new 
phase diagram for the SI transition which takes into ac- 
count the possibility of a coupling of the system to a 
dissipative bath. They argued that this coupling, which 
becomes important when the critical point is approached, 
can result in a new, metallic-like phase. In this picture, a 
direct SI transition is still possible for very weak coupling, 
while for a stronger coupling the system goes through a 
metallic phase and is truly superconducting only at the 
lowest magnetic fields. 

The fact that the typical sheet resistances of our sam- 
ples are about a factor of five higher than those in which 



G 



resistance leveling was observed |44j might just mean that 
our samples are in the weak coupling regime. However, 
the correlation length exponent determined in our exper- 
iment for the magnetic field-tuned transition, using two 
different methods, on different physical samples and at 
several levels of disorder was found to be v « 0.7, which 
is not consistent with the exponent expected from the 
classical 2D percolation theory, v = 4/3, even with much 
more generous error bars. 

A coherence length exponent of 0.7 is also inconsis- 
tent with what was believed to be an exact theorem, ]49| ] 
which predicts v > 1 in two dimensions in a presence 
of disorder. It is interesting to note that our exponent 
agrees with the result of the classical 3D XY model which 
is suggested to be relevant in the absence of disorder . 
Numerical simulations of a (2+l)-dimensional XY model 
H and the Boson-Hubbard model at T = (h]] without 
disorder also find z = 1 and v = 0.7. However, recently 
it was suggested that the nature of disorder averaging 
may introduce a new correlation length, different from 
the intrinsic one, which might lead to v < 1 even for a 
disordered system |6(| . 

There is also a possibility that the local dissipation cou- 
pled to the phase of the superconducting order param- 
eter due to gapless electronic excitations might change 
the universality class of the system and lead to a non- 
universal critical resistance [0. The critical resistance 
would then be expected to increase with increasing damp- 
ing due to dissipation. The latter would be expected to 
increase with decreasing normal state resistance. How- 
ever, we observe that the critical resistance decreases as 
the normal state resistance decreases, which is exactly 
the opposite of the behavior predicted by Wagenblast et 
al. 

We should note that v « 0.7 was also found for the 
magnetic field-tuned insulator-conductor transition in Si 
MOSFET samples |6l|| , suggesting a possible connection 
between the two phenomena. 

Our results for the magnetic field-tuned SI transition 
seem to be consistent with the predictions of bosonic 
models, rather than percolation models. This is fur- 
ther supported by the transport studies in the insulat- 
ing regime, where the magnetoresistance cannot be ex- 
plained by the weak localization theory only JO]] , and the 
temperature dependence of the resistance fits the predic- 
tions of Das and Doniach |6^] for the bosonic conduction. 
These observations still need to be reconciled with the 
results of the tunneling experiments which find no super- 
conducting gap in the insulating regime. The tunneling 
experiments might however be emphasizing regions of the 
samples containing quasi-localized single electron states 
below the gap, or those in which the amplitude fluctu- 
ations break the system into superconducting "islands" 
with finite spectral gaps in the density of states, as re- 
cently predicted |2C| ] . A highly non- uniform gap has also 
been predicted by Hcrbut |63| for the case of large dis- 
order. This problem might be clarified using spatially 
resolved scanning tunneling spectroscopy at low temper- 



atures, which may be able to detect local variations in 
the density of states. 

Such studies might also help answer the question as 
to why v different for the thickness- and magnetic field- 
tuned transitions on the same samples. In the case of 
the thickness-tuned transition, the correlation length ex- 
ponent is close to what might be expected from the per- 
colation theory. There is a major difference between the 
magnetic field-tuned and the thickness-tuned transitions: 
when the transition is tuned by the magnetic field, the 
microstructure of the sample stays fixed, while in the 
case of the thickness-tuned transitions it changes slightly 
with each film in the sequence. It may be that in this case 
the percolation effects become relevant, complicating the 
determination of the critical exponents fj4| . 

Finally, the shape of the phase boundary poses a fur- 
ther challenge to theorists. We are currently investigat- 
ing the role of the dissipation in this system in more 
detail, using a 2D electron g substrate, similar to 

the experiment of Rimberg et al. [j5(| . 
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FIG. 1. Resistance per square as a function of temperature 
in different magnetic fields, ranging from OkG (bottom) to 
12kG (top), with lkG increments. 
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FIG. 11. The critical resistance as a function of tempera- 
ture for a 12.353 A thick film. 

command, command. 

FIG. 3. Normalized resistance per square as a function of 
the scaling variable, T^ X ^ VZ \B — B c \. Each symbol represents 
one film at different temperatures (only a small portion of the 
data is shown for clarity). Inset: The fitting a power law to 
the temperature dependence of the parameter t determines 
the value of vz. 

FIG. 4. The critical exponent product vz for the magnetic 
field-tuned transition as determined by the inverse slope of 
dR/dB at the critical value of B c plotted vs 1/T. 

FIG. 5. Resistance per square as a function of magnetic 
field at different electric fields across the film: 0.5 (botttom), 
1.0, 1.5, 2.0, 2.5, 3.0 and 3.5 V/m (top). Only B < B c is 
shown where the resistance increases with increasing electric 
field. The temperature is 0.7K. 



FIG. 6. Resistance per square as a function of the scaling 
variable, t\B — B c \, for different electric fields: 0.5 , 1.0, 1.5, 
2.0, 2.5, 3.0 and 3.5 V/m. Here t = E~ 1/v(z+1) is treated as an 
adjustable parameter to obtain the best collapse of the data. 
Inset: The fitting a power law to the temperature dependence 
of the parameter t determines the value of vz. 

FIG. 7. Normalized resistance per square as a function of 
the scaling variable, t\d — d c \, in different magnetic fields: 
0.5 (squares), 1.0 (circles), 3.0 (crosses), 4.5 (triangles) and 
7.0 kG (diamonds). Inset: The fitting a power law to the 
temperature dependence of the parameter t determines the 
value of vz. 

FIG. 8. The critical exponent product vz for the thick- 
ness-tuned transition as determined by the inverse slope of 
dR/dd at the critical value of d c plotted vs 1/T. 

FIG. 9. The phase diagram in the d-B plane in the T=0 
limit. The points on the phase boundary were obtained 
from thickness tuned transitions (triangles) and magnetic 
field-tuned transitions (circles). The solid line is a power law 
fit. Here d c is taken to be the critical thickness in zero field. 

FIG. 10. The critical resistance as a function of the critical 
field for a series of bismuth films. Here R c decreases with 
increasing thickness, as thicker films have lower normal state 
resistances and higher critical fields. 



FIG. 2. Resistance per square as a function of magnetic 
field for a bismuth film close to the transition. Different curves 
represent different temperatures: 0.15, 0.17, 0.19, 0.2, 0.25, 
0.3 and 0.35K. 
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